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In quantum physics, all measured observables are subject to statistical uncertainties, which arise 
from the quantum nature as well as the experimental technique. We consider the statistical uncer- 
tainty of the so-called sampling method, in which one estimates the expectation value of a given 
observable by empirical means of suitable pattern functions. We show that if the observable can be 
written as a function of a single directly measurable operator, the variance of the estimate from the 
sampling method equals to the quantum mechanical one. In this sense, we say that the estimate 
is on the quantum mechanical level of uncertainty. In contrast, if the observable depends on non- 
commuting operators, e.g. different quadratures, the quantum mechanical level of uncertainty is not 
achieved. The impact of the results on quantum tomography is discussed, and different approaches 
to quantum tomographic measurements are compared. It is shown explicitly for the estimation of 
quasiprobabilities of a quantum state, that balanced homodyne tomography does not operate on 
the quantum mechanical level of uncertainty, while the unbalanced homodyne detection does. 
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I. INTRODUCTION 

The theory of measurements in quantum physics has 
a long tradition. First, it requires a good understand- 
ing of the measurement devices. In quantum optics, 
the foundations for the photodetector have been estab- 
lished in the pioneering works of Mandel [1] , Kelley and 
Kleiner [5j as well as Glauber (3J. On this basis, more 
elaborate measurement techniques have been developed, 
such as balanced homodyne detection [U [5] , unbalanced 
homodyne detection [B] and eight-port homodyne detec- 
tion schemes [7J [5] ■ All mentioned schemes provide the 
possibility to collect information to completely charac- 
terize an arbitrary quantum state of light. 

Once, the experimental techniques are well-known, one 
needs suitable tomographic methods to convert the raw 
experimental data into a convenient representation of the 
quantum state. In case of balanced homodyne measure- 
ments, in which one records phase dependent quadra- 
tures, there are several approaches which have been pro- 
posed. In [S], the so-called inverse Radon transform is 
applied to estimate the Wigner function of the quan- 
tum state. This approach can also be generalized to ob- 
tain different quasiprobability distributions, such as the 
Glauber-Sudarshan P function 10 1. In order to calcu- 
late density matrices, one may choose between Fourier 
techniques |llj . direct sampling schemes 12 13] or more 
involved maximum likelihood methods [151 116) . 

However, there has only been little interest in the sta- 
tistical properties of quantum state estimation. The the- 
oretical foundations have already been considered in |17j , 
but in particular for quantum state tomography, the 
statistical uncertainties of the estimates are not exam- 
ined deeply. First examinations have been performed 
in [TB] . The advantage of the sampling approach is that 
it directly gives an simple estimate of the statistical un- 
certainty of the estimated quantity [19], which we will 
shortly review below. In case of the maximum likelihood 
approach, statistical uncertainties can be evaluated by 



the inverse of the so-called Fisher information matrix, 
which requires some numerical effort [20j . Alternative 
methods for the estimation of uncertainties can also be 
found in 

In the present work, we compare the statistical un- 
certainty, which arises in the direct sampling approach, 
to the quantum mechanical variance, which provides a 
lower bound set by the quantum nature of light. In Sec. 
II, we consider functions of a single observable, which 
can be directly measured, and show that the variance of 
the corresponding pattern function equals to the quan- 
tum mechanical one. In Sec. Ill, we examine observ- 
ables for which complete quantum state tomography is 
required. Explicitly, we show that sampling from bal- 
anced homodyne detection data does not operate on the 
quantum mechanical level of uncertainty, and find that 
the statistical independence of quadratures at different 
phases is the reason. We also discuss the estimation 
of phase-space distributions in Sec. IV, and show that 
the unbalanced homodyne detection scheme provides es- 
timates with quantum mechanical uncertainties. Section 
V is dedicated to an example in quantum state tomogra- 
phy to illustrate the impact of the results. 

II. FUNCTIONS OF A SINGLE OBSERVABLE 

A. General considerations 

As a first step, let us assume that we observe a single 
physical quantity A and estimate the expectation value 
of a function of this observable, F(A). As F(A) is a 
Hermitian operator, it can be written in its spectral de- 
composition 

F(A)= I F(A)\A)(A\dA, (1) 

J A 

where \A) is an eigenvector of A with the eigenvalue A. 
The set A is the set of all eigenvalues. In case of dis- 
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crete eigenvalues, the integration has to be replaced by 
the corresponding sum. The eigenvectors can be chosen 
orthogonal, 



(A\A') = S(A-A'), 



(2) 



where the right hand side has to be interpreted as the 
Kronecker symbol in the case of discrete eigenvalues. Fur- 
thermore, the eigenstates provide a resolution of identity, 



\A)(A\dA = 1. 



(3) 



A 



Let us now consider an experiment which records a 
set of N eigenvalues {Aj}^ =l from A as outcomes. The 
underlying quantum state shall be denoted by p. Then, 
the empirical expectation value of some function F(A) 
can be estimated as 



1 N 



(4) 



The tilde indicates that this quantity is a random num- 
ber, since it is obtained from measured (and therefore 
random) values {Aj}^ =l , whose probabilities shall be de- 
noted by p(Aj). Therefore, the expectation value of this 
random number F is given by 



F= [ P (A)F(A)dA. 

J A 



(5) 



As above, the integral has to be seen as a sum over the 
probabilities if the set of eigenvalues is discrete. The key 
point is now that Eq. Q provides a good estimate for 
the quantum mechanical expectation value of the operator 
F(A). The probabilities of the outcomes are connected 
to the underlying quantum state by 

p(A)=Tr{p\A)(A\}. (6) 

Inserting this into Eq. ([5| and applying ([I]), we find 



F = Tr{pF(A)}. 



(7) 



The Eqs. Q,([5} and |7| form the basis of the sampling 
technique. In order to find an unbiased estimate of the 
expectation value of F(A), one simply has to insert his 
measurement outcomes {Aj}^ =1 into the so-called pat- 
tern function F(A) and calculate the empirical mean F 
of the values according to Eq. Q. 

For the ability to make justified statements, one 
still needs a measure of the uncertainty of the esti- 
mate Q. The empirical variance of the sampling points 
{/'■:. ' i is given by 



1 N 



Ff 



(8) 



This number quantifies the spreading of the points 
F{Aj). The factor j^zi guarantees that the estimate 
is unbiased, i.e. 



'F(A) 



p(A) F(A) - F 



dA. 



(9) 



Practically, the estimation of the empirical variance re- 
quires the calculation of the second moment of F(A), 



^ (2) = ^X>(^) 



(10) 



We easily see that this is exactly the sampling equa- 
tion Q, just with the square of the function F(A). 
Therefore, and due to the orthogonality of the eigenstates 
|j4), its expectation value again equals to the quantum 
mechanical one, 



F( 2 ) = Tr{pF (A)}. 



(11) 



The same holds for arbitrary moments of the function 
F(A). As a consequence, also the quantum mechanical 
variance of F (A), 

((AF(A)f) = Tr{p[F(i)] 2 } - [Tr{pF(i)}] 2 (12) 

can be estimated without bias from the set of sampling 
points {F(Aj)}jL 1 : 



a% (A) = ((AF(A)Y) 



(13) 



Finally, we are interested in the statistical uncertainty on 
the estimate Q. From estimation theory, this is just the 
empirical variance crf^) of the sampling points , divided 
by the number of measurements: 



2 _ 2 
°> — a F(A) 



IN. 



(14) 



The factor 1/iV guarantees that for an increasing num- 
ber of data points N, the uncertainty of the empirical 
mean value is decreasing. For N — ¥ oo, the latter ap- 
proaches stochastically the quantum mechanical expec- 
tation value. 

In conclusion, we may estimate the quantum me- 
chanical expectation value of the operator F(A) by the 
sampling equation Q, and the variance of F is ex- 
actly expected to match the quantum mechanical vari- 
ance {(AF(A)) 2 ), divided by the number of data points. 
In this sense, we may state that the determination of 
Tr{pF(A)} can be done on the quantum mechanical level 
of uncertainty. There are no other sources of noise con- 
tributing to the uncertainty, and it is not possible to 
achieve less fluctuations with classical statistical means. 

We emphasize that so far we only considered a function 
of the directly measurable operator A, making this result 
possible. After briefly discussing some examples, we show 
that the situation becomes completely different when we 
require several non-commuting observables to estimate 
the expectation value of an operator F. 
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B. Functions of a single quadrature 

As a first example, let us consider the measurement 
of the quadrature operator x, which is frequently done 
in balanced homodyne measurements. In this case, the 
set of eigenvalues is given by the continuous spectrum 
A = M. According to the above calculations, any func- 
tion of a single quadrature can be estimated at the quan- 
tum mechanical level of uncertainty. This includes all 
kinds of moments of the quadrature, for instance nor- 
mally ordered ones. However, note that we only may use 
a quadrature at a single phase. If we consider functions of 
quadratures of different phases, the situation will become 
completely different, as we will show below. 



C. Functions of the photon number 

Photon number resolving detectors can record the out- 
comes of the photon number operator n. Here, the set 
of eigenvalues is the discrete spectrum A = N. Again, 
we can state that the expectation value of any function 
of the photon number operator can be estimated at the 
quantum mechanical level of uncertainty. We still note 
that in practice one can only record a finite number of 
measurements, leading to some uncertainty on photon 
number probabilities with very low values, which typi- 
cally occur for very large photon numbers. However, this 
problem can be minimized by increasing the number of 
measurements . 



III. EXPECTATION VALUES FROM 
BALANCED HOMODYNE DETECTION 

Measurements of a single operator, such as the quadra- 
ture or the photon number, can not characterize a quan- 
tum state completely. For some operators, one needs 
more information about the state in order to estimate 
some expectation value. For instance, measurements of 
the quadrature distributions for all phases in [0, 7r) are 
informationally complete 23, 24], and we may calculate 
any expectation value from this measurement outcomes. 
On the other hand, it has been shown that so-called 
quasiprobability representations of quantum states can 
be retrieved by photon number resolved measurements, 
when one displaces the quantum state in phase space be- 
fore the measurement [6]. In the following, we consider 
these methods more in detail. 



operator p of the state and of the observable F. In gen- 
eral, the characteristic function of the Wigner function is 
defined as 



$ # (/3) = Tr{FD (/?)}, 



(15) 



with D(P) = e^ at 13 a being the well-known displace- 
ment operator [25] . For Hermitian operators, we have 
the relation $j=,(-/3) = $^(/3). Moreover, if F is the den- 
sity operator of a state, we will omit the index of <&(/?) 
throughout the paper. Conversely, the operator may be 
retrieved from its characteristic function by 



F = 



1 



d 2 ^U(3)D(f3). 



(16) 



With these quantities at hand, we may calculate the ex- 
pectation value of F with respect to the state p by 



Ti-{pF} = 



1 



d 2 /?$(^)$y/?). 



(17) 



Let us now express the second moment of F in terms of 
characteristic functions. Inserting Eq. ( 16 ) for the oper- 
ator F, we obtain 

Tr {pF 2 } = -L J d 2 ?' J d 2 0"^')^(P") 

xTr{pb{p')b{p")}. (18) 

By applying the equality 

D(a)b{P) = D(a + P)e am W\ (19) 

and writing both integrals in polar coordinates, we find 
the final expression 



Tr{pF 2 } 



1 



dip d<P db' db" \b' | \b"\Q(b'e iv + b"e lct >) 



7T 

x^* p {b / e zv )^* p (b /l e i4 ')e tb ' b " sin(,fi -' p >. (20) 



In our notation, the integrals over the angles <p, <f> range 
from to 7r, while the integration over b',b" covers the 
full real line. The variance can then be easily calculated 
as 



Tr{(AF) 2 } = Tr{pF 2 } - [Tr{^}] 2 



(21) 



Equation ( 20 ) will be the reference for comparison with 



the variance arising in sampling methods. 



A. Quantum mechanical expectation values 

For having a meaningful reference, we start with the 
calculation of the quantum mechanical variance of some 
operator F. Here, we express it in terms of the charac- 
teristic functions of the Wigner function of the density 



B. Balanced homodyne measurements 

1. Sampling from balanced homodyne measurement data 

Let us assume that the state p is sent to a balanced 
homodyne detector [5J [5], recording quadrature values 
{xj, (Pj)}j—i at different phases ipj € [0, it). Here, we 
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assume that the phase values are uniformly distributed 
within this interval. Then, the quadrature Xj follows the 
quadrature distribution p(xj;ipj), which is conditioned 
on the value of tpj. The joint probability distribution is 
given by p(x; ip) / w. 

Sampling is an established technique to estimate the 
expectation value of an Hermitian operator F from this 
set of data by an empirical mean of a suitable pattern 
function /(x, <p), 



Finally, the variance of the estimate F for the expectation 
value of the operator F can be obtained from 



/N 



(30) 



This number quantifies the uncertainty of the esti- 
mate J22|. 



Let us now examine if we can expect the empirical vari- 
ance Eq. ( 29 ) to match the quantum mechanical variance 



1 N 

3=1 



(22) 



Analogously to Eq. Q , this number F is a random vari- 
able, whose expectation value is given by 

F = [°° dx [" dv>?^f(x,<p). (23) 

J-oo JO K 

The pattern function has to be designed such that this 
expectation value equals to the quantum mechanical one, 



of the operator F. Contrarily to the procedure above, 
the operator F now depends on quadratures at different 
phases tp, whose operators do not commute. For getting a 
deeper understanding, we concentrate on the second mo - 
ment of f(x, tp). By using the inverse relation of Eq. (251, 



p(x;ip) 



1 

2^ 



we find 



2 _ 



dx 



dbe- lbx $(ibe^), 



dip [f{x,tp)\ 

7T 



F = Tr{pF}. 



(24) 



1 

2^2 



Let us now find such a suitable pattern function belong- 
ing to the operator F. The characteristic function of the 
state is connected to the quadrature distribution as 



dxdipdbdb' db" \b'\\b"\e l{b ' 



x^ibe iv )$%(ib'e iv )$%(ib''e i,p ). 



(31) 
(32) 

hb"~b)x 

(33) 



d>(/3) 



dxp (x]arg(fi)- §) e i|/3 K (25) 



Inserting this equation into Eq. (17) and writing the in- 
. over /3 in j 

Tr{pF} = I dx I dtp- 



tegration over /? in polar coordinates, we obtain 

p(x;tp) 



We substitute ip — ¥ ip — 7r/2 in order to remove the imag- 
inary unit i in the arguments of the characteristic func- 
tion. A careful analysis shows that we do not have to 
change the integration domain due to the periodicity of 
the integrand. Moreover, the integral over x can be eval- 
uated as 



oo Jo 

'CO 



1 

27 



dx e i ( b '+ b "-V x = 5(b' + b" - 6), (34) 



db \b\ e ibx <&*~ (ibe iLp ) (26) which can be used to evaluate another integral in (33): 

-oo 



From this relation, we easily see that the pattern function 
is given by 

/oo 
db\b\e ibx ^%(ibe^). (27) 
-co 

By construction, the expectation value of this pattern 
function with respect to the joint quadrature distribu- 
tions p(x] tp)/ir always equals to the quantum mechanical 
expectation value, see Eq. pi] ). 

Moreover, we are interested in the empirical variance 
of the single data points, which can be estimated from 
the experimental data by 



f(x,p>) 2 = - dipdb' db" \b'\\b"m(b' + b")e i<p ) 



x<S>*Jb'e ilp )<f>*Jb"e ltp ). 



(35) 



Together with Eq. ( 29 ) , we find the theoretically expected 



variance of the sampling method. 

We stress that this equation is not directly evaluated 
in practice, since the underlying quadrature distribution 
is unknown, and we only have a sample of quadrature 
measurements {(xj, ipj)}f = -\ ■ I nstead, we use the empir- 



ical variance given in Eq. (28). However, the theoreti- 



cal expectation is given by Eq. ( 29 ) in combination with 



Eq. ( 23 ) and ( 35 ) and the basis for all following consid- 



'/(*.¥>) 



N - 



1 N 



(28) 



3=1 



being completely analogous to Eq. ([8| . Consequently, the 
expectation value of this variance is given by 



(29) 



erations. 



2. Comparison with the quantum mechanical variance 

Let us compare the quantum mechanical vari- 
ance (21) with the variance expected from the sampling 
method ( 28 ) . Since the first moments of the sampling 
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method and quantum mechanics are equal by construc- 
tion, it is sufficient to examine the second moments 
of f(x; ip) and F. The expressions from the sampling 
method ( |35| ) and the quantum mechanical calculation 
(20) look quite similar, but they are different: in ( [35] ) , 
one integration over (f> is missing. A closer look reveals 
that if one replaces 



1 



dej) 



d<f>5(cp — 4>) 



(36) 



in the quantum mechanical expectation ( 20 1 , one finds 
the expression for the pattern function (35). Obviously, 



since we only observe the quadrature distribution p{x; ip) 
at a fixed phase ip, and all measured quadratures for 
different ip are stochastically independent, the "correla- 
tions" of the quadrature distributions between different 
phases ip and </> are not taken into consideration. This 
is due to the fact that p(x; ip) is not the joint distri- 
bution of all quadratures, whose marginals are the ob- 
served quadrature distributions. The definition of the 
joint distribution suffers from the problem of the non- 
commutativity of the corresponding quadrature opera- 
tors and is closely related to the different phase-space 
distributions. As a consequence, the quantum mechani- 
cal variance is not equal to the empirical variance of the 
pattern function. In this sense, the statistics of the bal- 
anced homodyne measurements is not at the quantum 
mechanical level, since the expectation value of F can 
not be estimated with quantum mechanical uncertainty. 

The definition of a joint probability of two quadratures 
suffers from the problem of the non-commutativity of two 
quadratures for different phases. In this context there 
appear similar problems as the well known ambiguity of 
the definition of phase-space distributions. 

We also note that the statement does not change when 
the examined operator is phase-insensitive, i.e. 



(37) 



In this case, it is sufficient to record the completely phase- 
diffused quadrature distribution 



1 



P P d{x) 



■p(x;<p)dip, 



(38) 



which can be easily done by choosing a uniform phase 
distribution of the local oscillator coherent state. This 
scheme may be practically more simple, since one does 
not have to record the phase values. However, it is still 
necessary to guarantee the uniform phase distribution 
in (38). Moreover, it does not bring any statistical bene- 



fit, since the variance ( 35 ) of the pattern function remains 
the same. Therefore, the same number of data points is 
required to obtain the same statistical precision. 

Finally, we might ask if it is possible to experimentally 
construct a bipartite state described by the characteris- 
tic function <&(/3' + /3"), which appears in Eq. p0|. The 



arguments /?' and /3" are assigned to each of the two sub- 
systems. If this was possible, one could use it for esti- 



mating the quantum mechanical variance from Eq. ( 20 1 



Practically, we would need to measure joint quadrature 
distributions 



p(xi,x 2 ; <fii,<P2) 



1 



(2tt)2 j_ 



db' 



db"e^ ( - ib ' Xl+b " X2 " 1 



x$(ib'e ivi +ib"e iV2 ), 



(39) 



which seems to require joint balanced homodyne mea- 
surements on the bipartite state. 

The problem is just that $(/?' + j3") does not refer to 
a physical state. This becomes clear when we examine 
the covariance matrix. The required moments can be 
obtained by taking the derivatives of $(/3' + (3") at /?' = 
ft" = 0. We note that all moments of the bipartite state 
can be expressed in moments of the state $(/3), since 



ptk fil fim f\n 

— $(8' + B")\ 



Qfik+m 



l+n 



(40) 



Therefore, if the state described by $(/?) has a quadra- 
ture covariance matrix 



V x a 



the bipartite state can be characterized by the matrix 

In order to describe a physical state, this matrix has to 
satisfy the nonnegativity condition [2 6) 



where 



= 



J 
J 



c 2 + in> o, 



and J = 



1 

-1 



(43) 



(44) 



However, the 3x3 minors of C'2 are always negative, e.g. 



C 



V x 

xp 
Vx 



C X p + i Vx 



v P 



c 



c 



Xp 



Xp 

v, 



-v x . 



(45) 



Therefore, Eq. (43) is violated by the covariance matrix 



C2 of the bipartite characteristic function $(/3'+/3"), and 
3>(/3' + P") can never correspond to a physical quantum 
state, which could be generated in an experiment. As 
a consequence, it seems unfeasible to estimate the vari- 



ance ( 20 ) on the quantum mechanical level. 



3. Coherent displacement of states and operators. 

Next, we consider a more general setting. First, we 
assume that we look at a family of observables ^(7), 



() 



(b„( R) 


/ r (f>[ R- rv +\ 




\ l -a) 

1 r 1 



FIG. 1: Scheme for overlapping a signal state 3>o(/3) with a 
coherent state. 



being constructed by a coherent displacement of some 
initial operator F: 

F( 7 ) = D(j)FD(- 7 ). (46) 

The characteristic function of these operators is given by 

%)M = ^(«e^- r l (47) 

Furthermore, we also apply the displacement to the state 
p, but we look at a suitable experimental realization. Let 
us send the initial state, described by &o((3), to a beam- 
splitter with real transmissivity t and reflectivity r, sat- 
isfying i 2 + r 2 = 1 (see Fig. fl|. At the second input 
of the beamsplitter, we irradiate a coherent state with 
amplitude ta/r, described by its characteristic function 



$coh(/3) = e^V^-lta/r^-g-ipr/a. ( 48 ) 
The resulting output state has the characteristic function 

$08; a, t) = $ (i/3)$coh(r-/3) (49) 
= ^{tp)e ta ' p - ta ^e- {1 - t2 ^ 2 l 2 . (50) 



We have chosen the amplitude of the coherent state such 
that the reflectivity does not appear in all calculations. 
Note that the transmissivity t can also be used for taking 
the detector quantum efficiency into account. It is well 
known that an imperfect detector with efficiency r/ < 1 
can be modelled by first mixing the input state with a 
fraction of 1 — r\ of vacuum and subsequently performing 
the measurement with an ideal detector. The transmis- 
sivity of the corresponding beamsplitter is given by ^Jr\. 
By applying Eq. (50 1 on the state <&(/?; a, i) with trans- 



missivity t' ~ y/fj and coherent state amplitude a' = 0, 
we find 



= $(/3; a, V^)- (51) 



Therefore, an imperfect detector can be simply taken into 
account by replacing the transmissivity t of the beam 
splitter by the effective transmissivity t^Jrj. In conse- 
quence, we may consider only ideal detectors. 



Of course, the expectation value of the pattern func- 
tion corresponding to F(a) with respect to the new out- 
put state is exactly the quantum mechanical expectation 
again. Let us look at the second moment of the pattern 
function, by inserting Eq. ( 47 ) and ( 50 ) into ( 35 ) : 



f(x^) 2 



1 

7T .;„ 



dip db 1 db"\b'\\b"\<b {t{b' + b"y*)e~ {l - t2){h ' +b ' ,)2/2 

J —OO J — OO 



(52) 



We focus on the following aspect: We want to estimate 
the quantum mechanical expectation value of the opera- 
tor -FX7) with respect to the initial state. For this pur- 
pose, we have two possibilities: 

1. We do not displace the initial state and omit the 
beamsplitter. Mathematically, this is given by 
t = 1 and a — 0. The expectation value of F(j) 
is obtained by choosing the pattern function corre- 
sponding to F(j), i.e. by suitable calculations after 
the balanced homodyne measurement. 

2. We displace the state by a — — 7/t. The calcula- 
tions after the measurement only require the pat- 
tern function for F, i.e. we set 7 = in Eq. (52 1. 



The expectation value of the sampling procedure is in 
both cases the same, namely the quantum mechanical ex- 
pectation. However, the uncertainties differ: In the sec- 
ond case, the initial state enters as &o(tf3), which is the 
state after exposition to losses with r\ = t 2 . Therefore, 
we may expect a worse result than in the first scheme, 
which only depends on the perfect state. This finding 
also holds when we consider imperfect detection: If the 
quantum efficiency of the balanced homodyne detector 
equals to we had to replace t — > ^Jfj2 in the first 
case and t — > t^/rjd in the second case. Therefore, the 
beamsplitter which displaces the initial state introduces 
unavoidable losses. This has consequences for the discus- 
sion of quantum state tomography. 
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IV. RELATION TO QUANTUM STATE 
TOMOGRAPHY 

In the case of quantum state tomography, one is inter- 
ested to find a complete representation of the quantum 
state. Here, we discuss the representation with means 
of quasiprobabilities. In many cases, these quasiproba- 
bilitics can be represented as the expectation value of a 
displaced operator F, which is phase-insensitive: 



P(a) = Tr{pD(a)FD(-a)}. 



(53) 



The coefficients F n in the Fock basis expansion of F, 



n=0 



(54) 



determine a specific quasiprobability. For instance, the 
Wigner function is obtained by choosing F n = 2/ir(— l) n , 
while the Q function arises from F n = 5 n fi/ir. More gen- 
erally, the coefficients of fi-ordered quasiprobabilities [57] 
can be written as 



dbbVL w (b)L n (b 2 



(55) 



with L n (x) being the nth Laguerre-polynomial. For the 
family of s-parameterized quasiprobabilities, one chooses 
f2(&;s) = e _ ' 1_s ) b / 2 , for nonclassicality quasiprobabili- 
ties one uses a suitable nonclassicality filter [28] [29] . 

We have discussed different techniques for estimating 
expectation values of the form (53 1. First, we notice 



that the so-called cascaded balanced homodyning tech- 
nique [30] does not provide any advantage over the stan- 
dard balanced homodyne detection. The former corre- 
sponds to the method 2 described in the previous section, 
where one displaces the state and estimates F for 7 = 0, 
while the latter is realized by the method 1. Obviously, 
the former suffers a reduction of the quantum efficiency 
by the transmissivity of the first beamsplitter. This is 
not affected by the fact that F is phase-independent, and 
recording of the phase values is not required in the cas- 
caded measurement. One can only improve the situation 
by choosing a beamsplitter with high transmissivity. 

Second, we note that balanced homodyne measure- 
ments followed by sampling of quasiprobabilities does not 
work on the quantum mechanical level of uncertainty. 
Therefore, we can not expect that this scheme is optimal 
for this purpose. Indeed, if one is interested in quasiprob- 
abilities, there is a better alternative: The unbalanced 
homodyne detection technique is based on a different in- 



terpretation of Eq. (53): First, the quantum state is dis- 
placed in the same way described in Sec. |IIIB 3[ After- 
wards, the expectation value of the phase-independent 
operator F is sampled from photon number measure- 
ments. Since this only requires the measurement of a 
single observable, namely the photon number, the esti- 
mation of the quasiprobability at a specific point a is per- 
formed on the quantum mechanical level of uncertainty. 



Provided that the balanced homodyne detector and the 
photon number resolved detector had the same quantum 
efficiency, we recommend to choose the latter one, since 
we expect it to give better results. 

Finally, we emphasize that the unbalanced scheme is 
optimal for the estimation of the phase-space represen- 
tation at fixed points a, but does not cover correlations 
between different points a, a'. Therefore, if one wants 
to estimate quantities which require the knowledge of 
the quasiprobability at different a, one can not expect 
to achieve this at the quantum mechanical level of un- 
certainty as well. For instance, quadrature distributions 
are better measured in balanced homodyne detection 
schemes. In this sense, the unbalanced measurement is 
optimal for the estimation of quasiprobabilities, but not 
always the best in different cases. 



V. EXAMPLE 

To demonstrate the difference of the statistical un- 
certainty in balanced homodyne and photon-number re- 
solved detection, let us consider the determination of a 
nonclassicality quasiprobability of a squeezed state |28l 
151] . Its variances are V x = 0.5 and V p = 2.0, whereas the 
variance of the vacuum state shall be V vac = 1. We will 
apply 

Sl w (P) = I 'u(p)u3{p + f3/w)d 2 (3 



(56) 



as a filter, with w(/3) = (2/7r) 3 / 4 e^l 4 . The prefactor 
guarantees that £l w (0) — 1, and the filter width is fixed 
with w = 1.8. Moreover, the pattern function (27) is 
defined by choosing <J> # (/3) = ^VL^e^ 2 / 2 e a " p ~ afi \ 

For the balanced homodyne detection scheme, we gen- 
erate a set of N = 100000 data points, each consisting 
of a pair (ifj,Xj). The phase values ipj are uniformly 
distributed in the interval [0, tt), whereas xj follows the 
quadrature distribution p(x; <pj) conditioned on the value 
of ifj. From this simulated set of data, we sample the 
nonclassicality quasiprobability together with its statis- 
tical uncertainty. For reasons of simplicity, we assume to 
have an ideal detector, i.e. 77 = 1. 

In case of the unbalanced homodyne detection scheme, 
we calculate the photon-number distribution p n together 
with its variance theoretically, the maximum photon 
number is restricted to n = 20. Then we derive the 
statistical uncertainty from the result by means of linear 
error propagation. 

Figure [2] shows the nonclassicality quasiprobability. 
We observe negativities for some real a, being signa- 
tures of the nonclassicality of the squeezed state, the 
minimum is achieved at a — 0.6 with P(0.6) = —0.31. 
Figures [3] and |4] show the standard deviations crj,(a) and 
a u (a), which are obtained from balanced or unbalanced 
homodyne measurements respectively. They are calcu- 



lated from Eq. ( 29 ) and ( 21 ) , each divided by the number 
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FIG. 2: Nonclassicality quasiprobability of a squeezed state 
with V x = 0.5 and V p = 2.0. 




Im (a) Re(ct) 

FIG. 3: Standard deviation from balanced homodyne mea- 
surements of the nonclassicality quasiprobability in Fig. [2] 

of measurements N. Obviously, they show a completely 
different behavior. The uncertainty from balanced homo- 
dyne detection is more than a factor of 3 larger than the 
one from the unbalanced technique, the exact difference 
depends on the point in phase space and on the examined 
state. In particular, at a = 0.6, the homodyne measure- 
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FIG. 4: Standard deviation from unbalanced homodyne mea- 
surements of the nonclassicality quasiprobability in Fig. [2] 
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FIG. 5: Comparison of the results for balanced and unbal- 
anced homodyne tomography. The blue shaded area corre- 
sponds to one standard deviation. 



ment provides <7f,(0.6) = 0.191, leading to an insufficient 
significance of the negativity of 1.6 standard deviations. 
Contrarily, we have cr„(0.6) = 0.010 in the unbalanced 
case, leading to a significance of about 31 standard de- 
viations. Therefore, in case of equal quantum efficiency, 
the unbalanced scheme proves to be much better than the 
balanced one. Fig. [5] illustrates this conclusion clearly. 



VI. CONCLUSION AND OUTLOOK 

We compared different approaches for estimating the 
expectation value of some physical quantity with respect 
to its statistical uncertainty. First, we showed that when- 
ever one can estimate a quantity of interest as a function 
of a single operator, which can be directly measured, then 
the estimate is on the quantum mechanical level of uncer- 
tainty, i.e. the empirical variance equals to the quantum- 
mechanical one. In practice, this works for quadratures 
and photon number measurements, for instance. 

However, for many operators a direct measurement is 
not known, and techniques for quantum tomography have 
to be applied. We considered sampling methods, which 
are applied to phase-dependent quadrature data from 
balanced homodyne measurements. We show that the 
cascaded balanced homodyne measurement has no sta- 
tistical advantage over the standard technique, although 
the first method only requires a phase-randomized lo- 
cal oscillator. Moreover, both methods do not operate 
on the quantum mechanical level of uncertainty. On 
the contrary, the unbalanced measurement technique can 
achieve the quantum mechanical level of uncertainty. 
Therefore, the latter shall be generally the best of the 
considered methods. 

We also identified the main reason for the difference 
between the variance from balanced homodyne measure- 
ment and the quantum mechanical variance: It is due to 
the fact that quadratures at different phases are always 
measured stochastically independently. Therefore, there 
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seems to be a lack of "phase correlations" in the exper- 
imental data. The severe question is how this problem 
affects other quantum state reconstruction methods, like 
maximum likelihood techniques |15l 116). In particular, 
it is unclear if the latter method is able to perform the 
estimation on the quantum mechanical level. 

Our results have direct implications to the different 
approaches for reconstructing quasiprobabilities or den- 
sity matrices of states. We showed the advantage of the 
unbalanced measurement with the example of a nonclas- 
sicality quasiprobability of a squeezed state. Therefore, if 
photon-counting devices with quantum efficiencies com- 



parable to balanced homodyne detectors will be available 
in the future, our findings suggest to prefer the unbal- 
anced scheme for quantum state estimation. 
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